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Abstract 

We make a general study of possibility of generating solar scale A Q and the CHOOZ angle 
U e s radiatively by assuming that they are zero at some high scale. The most general neutrino 
mass matrix leading to this result is determined in a CP conserving theory This matrix 
contains four independent parameters which can be fixed in terms of physical observables. 
The standard weak radiative corrections then lead to non-zero A and U e 3 without drasti- 
cally altering the other tree level results. As a consequence, both A Q and U e z are predicted 
in terms of other physically observable parameters. These predictions are insensitive to spe- 
cific form of the neutrino mass matrix. The solar scale and U e ^ are strongly correlated with 
the effective neutrino mass m ee probed in neutrinoless double beta decay. In particular, the 
LMA solution to the solar neutrino problem arise for m ee close to the present experimen- 
tal limit. An example of specific texture is presented which predicts maximal atmospheric 
mixing and tan 2 9 & ~ 0.5 for the solar mixing angle 6 & . 
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Introduction: Neutrino oscillation experiments have provided significant information 
on neutrino masses and mixing [JTJ] . The " standard" picture emerging from analysis of various 
experiments is the existence of two hierarchical (mass) 2 differences and two large and one 
small or zero mixing angle among three neutrinos. The overall scale of neutrino masses 
is not fixed directly by neutrino oscillation experiments. This complimentary information 
can be obtained from direct neutrino mass measurements and from neutrinoless double 
beta decay (0vf3j3) experiments ||. These types of experiments have so far provided only 
upper limits. The neutrino mass (assuming no mixing) is constrained by tritium f3 decay 
experiments to be less than 2.2 eV. The best limit from 0uj3(3 decay experiments is 



X^ei^J < 0.38 heV at 95% CL , (1) 



where h ~ 0.6 — 2.8 denotes the uncertainty in nuclear matrix element |Q. U denotes here 
the neutrino mixing matrix and m Vi [i = 1,2,3) are neutrino mass values which can take 
either sign. 

While the hierarchical neutrino masses cannot be ruled out at present, the presence of 
large mixing angles hints at an almost degenerate pair of neutrinos. This will become a 
necessity if m ee would be found to be significantly larger than the atmospheric scale. This 
would require all three neutrinos to be nearly degenerate || if mixing among them is also 
to account for the solar and atmospheric neutrino results. The mass patterns with two || 
or all the three |5],[7|,[8| nearly degenerate neutrinos are therefore of considerable importance. 

If overall scale of neutrino masses is larger than the atmospheric scale then one would like 
to understand why neutrino (mass) 2 differences (particularly, the solar scale A ) are much 
smaller? Interesting possibility is to assume degenerate neutrinos which get split by 
radiative corrections PJIOf induced through charged current interactions in the standard 
model (SM) or in the minimal supersymmetric standard model (MSSM). Advantage of 
this scheme is its predictive power. The most general 3x3 neutrino mass matrix with 



completely degenerate spectrum is characterized JO]] in terms of two angles and one phase. 
These three parameters determine all neutrino masses and (complex) mixing after the known 
radiative corrections are included. Unfortunately, this predictive scenario does not give M 
phenomenologically required description of neutrino masses and mixing. 

The next best possibility is to assume that only two of the three neutrinos are exactly 
degenerate at high scale. The radiative corrections then lead to the solar splitting within 
this scheme. It is possible to do a quite general and fairly model independent analysis of this 
situation which we present in this paper. We assume that neutrino masses, the atmospheric 
scale and two large mixing angles are tree level effects and are already described by neutrino 
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mass matrix specified at a high scale. We require this mass matrix to have vanishing solar 
scale and vanishing CHOOZ |JTI| angle in flavour basis. Neutrino mass matrix with such 
property can be characterized by four independent parameters in a CP conserving theory. 
Solar scale and the CHOOZ angle U e 3 arise after radiative corrections and represent generic 
prediction of this scheme. These predictions are found to be model independent and hold 
for all the matrices under consideration. We now present this analysis and discuss its 
consequences. 

General analysis of pseudo Dirac neutrinos: Let us consider a CP conserving theory 
specified by a general 3x3 real symmetric neutrino mass matrix M u0 . This matrix can 
always be specified in flavor basis corresponding to diagonal charged lepton masses. We 
adopt the following general parameterization for M u q in the flavor basis. 



M, 



vO 



( \ 

Si t u 
t s 2 v 



(2) 



yu v s 3 j 

We assume that the above M u0 describes physics at a high scale M X - M u0 is required to 
yield vanishing solar scale and CHOOZ angle at Mx- Let us derive conditions on elements 
of M u0 for this to happen. 

Vanishing of the solar scale requires that two of the eigenvalues of M u0 are degenerate 
with masses (m, m) or (m, — m) (m > 0). The relative angle between the degenerate pair 
can be rotated away in the former case in a CP conserving theory. This is not true in case 
of masses differing in their sign. Thus barring possibility of radiative amplification Jl3| , the 



former case will not lead to large solar angle and we concentrate on the second possibility 
with masses (m, — m). Such a pair is equivalent to a Dirac neutrino invariant under some 
global U(l) symmetry. The standard weak current would violate this symmetry in general 
nj and the Dirac state would split into a pair of pseudo Dirac neutrinos. General conditions 



under which this happens in case of three generations were discussed in [F5]. In particular, 
the M u0 should satisfy 

tr(M M )J2\ = detM M , (3) 

i 

where Aj represents the determinant of the 2x2 block of M u q obtained by blocking the i f h 
row and column. 

When the above condition is satisfied, eigenvalues of M u0 are given by (m, —m, T) where 



m =J-J2 A i T = tr(M u0 ) 



Let Uq diagonalize M, 



Uq M u q Uq = Diag.(m, —m, T) 



(4) 



Since M u0 is specified in the flavor basis, Uq defined above represents physical neutrino 
mixing matrix at tree level. The electron neutrino survival probability in reactor experiments 
such as CHOOZ is given by (£/o) e 3 which we require to be zero. One can show that M v q 
satisfies eq.(|3|) and leads to (?7o)e3 = provided 



t 



Oi + s 2 )(si + a 3 ) ; 
uv 



(5) 



S\ + s 2 

The second equation does not hold in a special case with v — 0. In this case t and u are 
unrelated and the above two conditions uniquely lead to the following M u0 : 



M, 



;/0 



/ \ 

s t u 
t -s 



i u — s j 



(6) 



The detailed phenomenological consequences of the special solution given in eq.@ were 
worked out in [0. Here we consider rest of the the neutrino mass matrices specified by 
restriction given in eq.(|5|). Let us parameterize the mixing matrix Uq by 



Uq 



c 0o s 4>o 



\ 



(7) 



The mixing angles are determined using eq.(||) and eq.(|): 



tan i 



\m — s\ 
m + Si 



tan 6*n = — . 

t 

Here, m denotes the common mass of the Dirac pair and is given by 



(8) 



m = <J s\ + t 2 + u 2 



We note that 
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• The effective neutrino mass probed in the Qp(5f3 is given by 

m° ee = st 

• At the tree level, there is only one (mass) 2 difference which provides the atmospheric 
scale 

A 0A = \m 2 -T 2 \. (9) 

Corresponding mixing angle (= 9 A ) coincides with 9q and is large when t ~ u: 

sin 2 29° A = sin 2 26 . (fO) 

• There is no solar splitting at this stage but would be solar mixing angle 9q coincides 
with (J) and is given by 

tan 2 ^ = ^^p (11) 
° m + m° ee 

The above relations are valid for the most general M M with parameters satisfying eq.(^). 
They are derived at tree level but as we demonstrate latter, radiative corrections do not 
substantially change them. The major effect of radiative corrections is to generate the solar 
splitting and a non-zero value for U e ^. It turns out that these quantities are not arbitrary 
but are predicted in terms of other observables irrespective of the detailed form of M v q. 
This happens because two conditions in eq.@ leave us with four independent parameters. 
They can be determined in terms of four observables namely, atmospheric scale and angle, 
effective mass probed in 0v(3(3 and the solar angle. The solar splitting and generated 
radiatively then no longer remain arbitrary but are determined in terms of these observables. 

We can express all parameters of M u0 in terms of observables using conditions of eq.(|): 

m =\^k%\ t 2 = cos 2 9° A (m 2 -m 02 ) 

T 2 = (m 2 - A A , m 2 + A A ) (for m 2 > A A , < A A ) (12) 



s 3 = cos 2 9° A T - sin 2 9° A m° ee s 2 = sin 2 9° A T - cos 2 9° A m° ee 



The solar splitting can be obtained (15| using the relevant renormalization group equa- 
tions @,fIIJ. The consequences of these equations have been discussed in a number of papers 

MM- 
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The radiatively corrected neutrino mass matrix is given by 



where 



M,. 



( i \ 
U 



11 



u 



M, 



it 

J-e 











1 

T 2 





v° 





1 

T 2 



(13) 



II = 1 + 5 a , 



with 



. . (14) 

•4nv' M z v ; 

M x here corresponds to a large scale and we take M x ~ 10 16 GeV; C — 2 5 cos^ j3 

respective 

cases of the standard model [] and the minimal supersymmetric standard model and 
a = e, /i, r. I g>t are calculable coefficient summarizing the effect of the gauge and the top 
quark corrections. 

Apart from the overall factor I g I t , the radiative corrections are largely dominated by the 
t Yukawa couplings and it is easy to determine neutrino mixing angle and masses keeping 
only 5 T corrections and working to the lowest order in 5 T . We now have 



U T M v U = Diag. (m^, m^, m„ 3 ) , 



with 



m,. 



I g I t (m + 5 T sin 2 6 A (m - m°J) + 0(5 2 T ) , 
I g I t (- m - 5 T sin 2 9° A (m + m°J) + 0(S 2 T ) , 
I g I t (T + 25 T T cos 2 8 A )+O(5 2 T ) , 

where we have used eq.(|l2]). The tree level mixing matrix Uq gets modified to a general U: 

( \ 



m V2 ^ 



m V3 



(15) 



U 



C^SqS^ + CgS^ CgC^ - S^SgS^ -SgC w 
\ -C^CgS^ + S^Sg S^S^Cg + SgC^ CgC^ J 



(16) 



1 The value of c in case of the standard model is given by 1/2(3/2) according to calculations in || 
( 0). We will use the value as in (Hj. 
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As before, the angles 0, 9 correspond respectively to solar and atmospheric mixing angles. 
These are now given by 

tan0 A = tan0^ (1 + / r _ 2 (m 2 + T 2 - 2Ts x ) ) +0(5 2 ), 
tan 2 6> =tan 2 ^°+0(5 2 ), (17) 

where #^ (eq. flToD ) and 9% (eq. (PD ) are tree level atmospheric and the solar mixing angles 
respectively. Note that the solar mixing angle does not receive radiative corrections to the 
lowest order in S T . 

The effective neutrino mass probed in Qv[3(3 is now given by 

m ee = I g Itm° ee = Igl t st. (18) 

The atmospheric scale also receive radiative corrections and is now given by 

A A = i(r< + m 2 J - m 2 3 = J 2 / 2 (A 0A + 25 T (m 2 sin 2 9 - 2T 2 cos 2 9)) + 0(5 2 T ) . (19) 

It is seen that all the tree level predictions receive small radiative corrections. Thus all 
the neutrino mass matrices leading to two degenerate states and characterized by eq. (|5|) are 
stable against radiative corrections. This is to be contrasted with the case of fully degenerate 
neutrino spectrum which leads to radiative instability in some specific cases 0. The non- 
trivial effect of the radiative corrections is generation of the solar splitting and a non-zero 
U e3 : 

A© = m - m w 4m ee d T — sin 9 A , 

cos 26 Q 



5 T T 'sin 2$ ' A \ /m 2 — m 2 

\u*\ = K\ ~ I ^ -| • (20) 

The above equations relate the solar scale and CHOOZ angle to other experimentally 
determined quantities as can be seen using eq.(|l"2D. Eq . (|20|) therefore represent basic predic- 
tion of the scheme defined by eq. (JH) . It is remarkable that all these matrices lead to unique 
predictions in eq.(^) which are insensitive to specific texture of the neutrino mass matrix. 

Let us now explore consequences of eq. ( pO"D . The atmospheric mixing angle and scale are 
experimentally well-determined: w (1.5 — 5) • 10 -3 eV 2 and sin 2 29 A ~ 0.8 — 1. The solar 



scale A and mixing 9 & are also highly constrained particularly after |18| the recent 
neutral current results from SNO ||19|| . Based on the global analysis of all the solar data, 
the only solutions allowed at 3a level are the large mixing angle solution (LMA) and the 



LOW solution with A ~ 10 7 eV 2 . The allowed ranges of parameters in these cases at 3a 
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are given approximately by by A Q « 3 ■ 10~ 4 - 2 ■ 10~ 5 eV 2 and tan 2 9 Q ~ 0.2 - 0.9 m case 
of the LMA solution and A Q « 3 • 10~ 8 - 1 ■ 10~ 7 eV 2 and tan 2 9 Q ~ 0.4 - 0.9 in case of the 
LOW solution. Both the small mixing angle and vacuum solutions are excluded at 3a. In 
particular, the solar mixing angle is found to be less than 45° in all the preferred solution a 
fact which plays an important role in the following. 

The predicted value of A Q and U e ^ are different in SM and MSSM due to different values 
of 5 T in these two cases. In case of SM, 5 T ~ 10 -5 while it can become larger for MSSM 
due to presence of tan/5. More importantly, sign of 5 T is different in these two cases. The 
negative values of 5 T in case of MSSM makes it unsuitable for the description of the solar 
data as we now argue. 

In analyzing the solar data, A is chosen positive by convention and mixing angle 9 Q is 
allowed to be greater than 45°. The sign of A as defined by eq.(|(J) is governed by the sign 
of m ee and 5 T . The sign of m ee also determines magnitude of the solar angle through eq.(|Tl"l). 
Positive (negative) values of m ee gives a 9 & less (greater) than one. Since S T is negative in 
case of the MSSM one needs negative m ee to obtain positive A Q with the result that tan 2 9 & 
becomes greater than one 0. Since the solar neutrino results do not allow tan 2 9 & > 1, 
MSSM radiative corrections as a mechanism to generate the solar splitting is disfavoured in 
the present context. In contrast, the SM radiative corrections give tan 2 9 Q < 1 as required 
for these solutions. We discuss this case now. 

The numerical value of A Q and U e % are correlated both with m ee as well as with the solar 
mixing angle. We show this correlation in Fig.(l) which displays variation in lO 5 -^- (solid) 
and 10 2 [/e3 (dotted) with m ee for typical values of tan 2 9 & needed for the LMA and LOW 
solutions. It is seen from the figure that LMA solution can be obtained for relatively larger 
value of m ee typically m ee > 0.1 eV. The minimum required value of m ee increases with 
decrease in tan 2 9 & . The LOW solution require much smaller but experimentally accessible 
values |20| of m ee around 0.05 eV. 

The predicted values of are generally smaller than the present limits as well as 
possible detection limit around ~ 0.05 for most ranges in the parameters. But if tan 2 9 & 



is ~ 0.5 — 0.8 then U e s is predicted to be in the range 0.01 — 0.1 and is correlated with the 
LMA solution. 



An interesting consequence |22fl of the correlation between LMA solution and large m e 



2 For positive m ee one needs to reverse the role of v\ and vi- The relevant tan 2 6q is inverse of 
eq . (|TT|) and is also greater than 1. 
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is as follows. For tan 2 6 Q ~ 0.2 — 0.7 and m e 
the degenerate pair lies in the range 



0.3 eV, the common mass m 



in,. 



I cos 20r. 



Of 



m « 0.5 - 1.8 eV . 

The third mass T is also required to be close to m since A^ ~ \m 2 — T 2 \. As a consequence, 
the LMA solution in these models is automatically correlated with almost degenerate neu- 
trino mass spectrum with a common mass much larger than the atmospheric scale. 

The above discussion is based on general class of matrices leading to pseudo-Dirac neu- 
trino. We now supplement this with a discussion of a specific texture. 
Almost degenerate neutrinos: Consider the following specific texture: 



M, 



i/0 



1 + e -2 
-2 1 - e 



2 
2 



1 



(21) 



The above texture is determined by only two parameters s and e. It satisfies conditions in 
eq.(||). Thus it leads to two degenerate eigenvalues and vanishing at a high scale. The 
eigenvalues of the above matrix are given by (m, —m, s(3 — e)) with 



m 



sV9 + 2e + e 2 



(22) 



It is seen that all neutrinos are degenerate when e = 0. It is known || that matrix with 
degenerate neutrinos cannot lead to the required mass pattern after radiative corrections. 
This is remedied here by introduction of a small e which leads to the atmospheric neutrino 
splitting at high scale: 



A. 



A 



AO 



8es 2 



(23) 



The above specific texture has four predictions. As in the general case, the radiatively 
generated solar scale and U e3 are predictions of the model. In addition, both the solar and 
atmospheric mixing angles instead of being arbitrary are fixed here by the specific texture. 

The solar splitting follows from the general expression in eq.(EH): 



A ~ 25 T s 2 (l + e) v / 9 + 2e + e 2 



(24) 



Two parameters s and e get determined by the values of A and A^. In particular, relatively 
large A ~ 10 _5 eV 2 needs small e and large s. The solar mixing angle is obtained using 
eqs. (0,P3) and is given in the small e limit by 
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tan 2 6> ^0.5 + 0(e) (25) 

Thus this texture automatically predicts large mixing angle which is in the range required 
for the LMA or LOW solutions. The atmospheric mixing angle is predicted to be maximal. 
Eq. ( pip therefore provides an example of the bi-large mixing patterns with almost degenerate 
neutrinos. 

One can determine required value of s, e from eqs.(|23],|24|). Choosing A ~ 5 • 10~ 5 eV 2 
and A^ = 3 • 10~ 3 , one finds 



e~ 1.8-10" d ; s = 0.45eV . 

Thus one can obtain A around 10 _5 eV 2 provided the effective neutrino mass m ee = s ~ 
0.4 eV. All the three neutrinos are almost degenerate with a common mass 1.3 eV which is 
not very far from the present experimental limit @. 

Summary: We discussed possibility of explaining small values of the solar scale and the 
angle U e ^ through radiative corrections by assuming that these are zero at a high scale. 
Restrictions to be satisfied by neutrino mass matrix for this purpose in the flavor basis were 
determined (eq|J). Since neutrino mass matrix can always be expressed in the flavor basis, 
eq.(|5|) provides general conditions for vanishing of U e 3 and solar scale in any model. 

We showed that the standard weak radiative corrections lead to the solar splitting re- 
quired on phenomenological grounds. Both A and U e % are predicted in terms of other 
observables. These predictions are remarkably independent of detailed form of M v q and 
remain true for any M u0 satisfying eq.(|). 

Detailed analysis presented here shows that one can obtain the most preferred LMA 
solution for m ee > 0.1 eV. Thus verification of LMA solution and moderate improvement 
in limit on m ee can rule out the entire class of solutions proposed here. The LMA solution 
also gets correlated in these models with almost degenerate neutrino mass spectrum and 
measurably large U e s ~ .01 — 0.1. 

It is not possible to obtain the correct solar parameters if radiative corrections are in- 
duced in MSSM. It may be possible to make MSSM also viable by allowing some nonzero 
£/e3 at high scale and/or by invoking additional sources of radiative corrections . 



Acknowledgments: I thank Saurabh Rindani for his interest in this work and for helpful 
discussions. 



10 



[1] For reviews of various models, see, A. Yu. Smirnov, hep-ph/9901208 ; R. N. Mohapatra, hep 



ph/9910365|; A. S. Joshipura, Pramana, 54 (2000) 119. 



[2] J. Bonn et al, Nucl. Phys. Proc. Suppl. 91 (2001) 273; Extensive discussion can be found in 
Y. Farzan, O. Peres and A. Yu. Smirnov, [hcp-ph /0105105| . 

[3] H. V. Klapdor-Kleingrothaus et al, Euro. Phys. Journal, A12 (2001) 147. 

[4] F. Feruglio, A. Strumia and F. Vissani, tiep-ph/0201291 . 



[5] D. Caldwell and R. N. Mohapatra, Phys. Rev. D48 (1993) 3259; A.S.Joshipura, Zeit. fur 
Physik, C64 (1994) 31; A. Ioannissyan and J. W. F. Valle, Phys. Rev. (1994); E. Ma, Phys. 



Rev. Lett. 83 (1999) 2514; E. J. Chun and S. K. Kang, Phys. Rev. D64 (2001) 09790 and [hep 
ph/9912524| ; E. Ma and G. Rajasekaran, Phys. Rev. D64 (2001) 113012; E. Ma, Phys. Rev. 
Lett., 83 (1999) 2514 and |hep-ph/0201225| ; R. Adhikari, E. Ma and G. Rajsekaran, Phys. Lett. 



B486 (2000) 134; M. Fukugita, M. Tanimoto and T. Yanagida, Phys. Rev. 57 (1998) 4429; M. 
Tanimoto,Phys. Lett. 483 B (2000) 417 and Phys. Rev. 59 (1999) 017304; H. Minakata and 
O. Yasuda, |hep-ph/9609276| ; F. Vissani, |hep-ph/9708"48^ ; E. J. Chun and S. Pokorski, Phys. 



Rev. D62 (2000) 053001. 

[6] R. Barbieri et al, Phys. Lett. B445 (1999) 239; A. S. Joshipura, Phys. Rev. D 60 (1999) 
053002; R. N. Mohapatra et al Phys. Lett. 474 (2000) 355; A. S. Joshipura and S. D. Rindani, 
Phys. Lett. B 464 (1999) 239 and Euro. Phys. Journal, C14 (2000) 85; K. S. Babu and R. N. 



Mohapatra, |hep-ph/020lT7e 



[7] J. Eliis et al, Eur. Phys. Jou. C9 (1999) 310; J. Ellis and S. Lola, Phys. Lett. 110 458 (1999) 
389. 

[8] A. Casas et al, Nucl. Phys. B556 (1999) 3; Nucl. Phys. B569 (1999) 82; Nucl. Phys. B573 
(2000) 652; Journ. of High Energy Physics, 9909 (1999) 015; A. Dighe and A.S. Joshipura, 
|hep-ph/0010079| ; N. Haba Y. Matsui and N. Okamura, Euro. Phys. Journ. C17 (2000) 513; 
N. Haba et al, Phys. Lett. B 489 (2000) 184; N. Haba et al Prog. Theor. Phys. 103 (2000) 
145; N. Haba and N. Okamura, Eurp. Phys. Journ. C14 (2000) 347; E. J. Chun, Phys. Lett. 



11 



B505 (2001) 155. 



[9] K. S. Babu, C. N. Leung and J. Pantaleone,Phys. Lett. B139 (1993) 191;P. H. Chankowaski, 
Z. Pluciennik,Phys. Lett. B316 (1993) 312. 

[10] S. Antusch et al, Phys. Lett. B519 (2001) 238 and Phys. Lett. B525 (2002) 130. 

[11] G. C. Branco, M. N. Rebelo and J. Silva-Marcos, Phys. Rev. Lett. 82 (1999) 683. 

[12] CHOOZ collaboration!, M. Apollonio et al, Phys. Lett. B466 (1999) 415. 

[13] K. R. S. Balaji et al Phys. Rev. Lett. 84 (2000) 5034 and Phys. Lett. B481 (2000) 33; Gautam 
Dutta jhlj>ph/02020971 and |hep-ph/0203222| ; N. Haba, N. Okumura and M. Sugiura, Prog. 
Theor. Phys. 103 (2000) 367. 



[14] L. Wolfenstein, Nucl. Phys. B 186 (1981) 147; S. T. Petcov, Phys. Lett. 110 B (1982) 245; 
S. T. Petcov and C. N. Leung, Phys. Lett. 125 B (1983) 461. 

[15] Anjan S. Joshipura and S. D. Rindani, Phys. Lett. B494 (2000)114. 



[16] A. S. Joshipura and S. D. Rindani, [hep-ph/ 0202064 



[17] G. L. Fogli et. al, |hep-ph/0106247| ; J. N. Bahcall et. al, |hep-ph/0106258| ; A. Bandyopadhyay 
et. al, hep-ph /0106264 ,P. Krastev and A. Yu. Smirnov, |hep-ph/0108177j ; V. Barger et. al, 



hep-ph/0106207 



[18] A. Bandyopadhyay et. al, |hep-ph/0020428| 6; V. Barger et al, |hep-ph/0204253| ; J. N. Bahcall 



et. al, |hep-ph/02043l4 



[19] SNO collaboration, Q. R. Ahmed et al |nucl-ex/0204008! and |nucl-ex/0204009 



[20] H. V. Klapdor-Kleingrothaus et al Modern Phys. Lett. A 16 (2001) 2409; C. E. Aalseth et 



al, |hep-ex/0202"0i~8t H. L. Harney, jhep-ph/0205293| and H. V. Klapdor-Kleingrothaus, |hcp" 



ph/020522§. 



[21] ICARUS collaboration, F. Arneodo, et al [hep-ex/010601£ 



[22] An incomplete list of the references is given by, H. V. Klapdor-Kleingrothaus, H. Pas and A. 



12 



Yu. Smirnov, Phys. Rev. D63 (2001) 73005; H. Minakata and O. Yasuda, |hep-ph/9609276| W. 



Rodejohann, Nucl. Phys. B597 (2001) 110 and |hep-ph/0203214| ; F. Vissani, |hep-ph/9904349| ; 



S. M. Bilenky et al, Phys. Lett. B465 (1999) 193; S. M. Bilenky, S. Pascoli and S. T. Petcov, 
Phys. Rev. D64 (2001) 053010; H. Minakata and H. Sugiyama, Phys. Lett. B532 (2002) 275 
and Phys. Lett. 526 (2002) 335; S. Pascoli and S. T. Petcov, frep-ph/0205"022] N. Haba and 
T. Suzuki, Mod. Phys. Lett. A17 (2000) 865; 

[23] P. H. Chankowski et al Phys. Rev. Lett. 86 (2001) 3486. 



13 



100 



10 : 




14 



